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ABSTRACT 



We measure the black hole mass in the nearby active galaxy Centaurus A (NGC 5128) using a new method based on spectroastrometry 
of a rotating gas disk. The spectroastrometric approach consists in measuring the photocenter position of emission lines for dilferent 
velocity channels. In a previous paper we focused on the basic methodology and the advantages of the spectroastrometric approach 
with a detailed set of simulations demonstrating the possibilities for black hole mass measurements going below the conventional spa- 
tial resolution. In this paper we apply the spectroastrometric method to multiple longslit and integral field near infrared spectroscopic 
observations of Centaurus A. We find that the application of the spectroastrometric method provides results perfectly consistent with 
the more complex classical method based on rotation curves: the measured BH mass is nearly independent of the observational setup 
and spatial resolution and the spectroastrometric method allows the gas dynamics to be probed down to spatial scales of ~ 0.02", 
i.e. l/IO of the spatial resolution and ~ 1/50 of BH sphere of influence radius. The best estimate for the BH mass based on kinematics 
of the ionised gas is then \og(MBH sinr/MQ) =^ 7.5 ± 0.1 which corresponds to Mbh = 9.6^^| x 10' Mg for an assumed disk incli- 
nation of i = 35°. The complementarity of this method with the classic rotation curve method will allow us to put constraints on the 
disk inclination which cannot be otherwise derived from spectroastrometry. With the application to Centaurus A, we have shown that 
spectroastrometry opens up the possibility of probing spatial scales smaller than the spatial resolution, extending the measured Mgu 
range to new domains which are currently not accessible: smaller BHs in the local universe and similar BHs in more distant galaxies. 
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1. Introduction 

One of the fundamental open questions of modern astrophysics 
is understanding the formation and evolution of the complex 
structures that characterize the present-day universe such as 
galaxies and clusters of galaxies. Understanding how galaxies 
formed and how they become the complex systems we observe 
today is therefore a major theoretical and observational effort. 

There is now strong evidence for the existence of 
a connection between supermassive black holes (hereafter 
BHs), nuclear activity and galaxy evolution revealing the 
so-called co-evolution of black holes and their host galax- 
ies. Such evidence is provided by the discovery of "relic" 
BHs in the center of most nearby galaxies, and by 
the tight scaling relations between BH masses {Mbh ~ 
10^ - IO'^Mq) and the structural parameters of the host 
spheroids like mass, luminosity and stellar velocity disper- 



Send offprint requests to: A. Gnerucci 



sion (e.g. Kormendv & Richstone 1995, Gebhardt et al. 200^ 
iFe rrarese & M erritt 200 0, .Marconi & Hunt,2003. .Haring & Ri 
2004 iFerrares e & Fordl l2005L lGrahamll2008l) . Moreover, while 
it has long been widely accepted that Active Galactic Nuclei 
(AGN) are powered by accretion of matter on a supermas- 
sive BH, it has recently been possible to show that BH 
growth is mostly due to accretion of matter during AGN ac- 
tivity, and therefore that m ost galaxies went through a phase 
of strong nuclear activity ( Soltan] Il982l lYu & Tremaind l2002l 
Marconi et al..,2004) . It is believed that the physical mechanism 
responsible for this coevolution of BHs an their host galax- 
ies is probably the feedback by the AGN, i.e. the accreting 
BH, o n the host galaxy (Silk & Rees 1998, Fabian & Iwasa;^ 
1999, Granato et alj|2004j, IDi Matteo etal J 120051 IMencil 120(361 
Bower et al...2006) . 

The clearest sign of co-evolution, the scaling relations be- 
tween BH masses and host galaxy properties should be then se- 
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cured by increasing number, accuracy and mass range of existing 
measurements. 

Supermassive BHs are detected and their masses measured 
by studying the kinematics of gas or stars in galaxy nuclei and, 
currently, there are about ~ 50 BH mass me asurements mos t 
of which in the ~ 10^ - lO^M© range (e.g. Sani et al.ll20"T()h . 
The majority of these measurements are made with longslit 
spectroscopy, but the development in recent years of Integral 
Field Unit (hereafter IFU) spectrographs has allowed some im- 
provements. IFUs have proven to be powerful tools to study 
galaxy dynamics as they provide two dimensional coverage of 
the source without the restrictions of longslit spectrographs, 
also plagued by unavoidable light losses. Recent studies have 
presented measurements of BH masses in galactic nuclei us- 
ing i ntegral field spectroscopy of gas o r stell ar spectral fea tures 
(e.g. iDavies e t al.''20Q6'. 'Nowak et al."2007', 'Nowak et al. 2010' 
Krainovic et al.i i2007. .Krainovic et al., ,2009, Cappellari et al.i 



120091 lNeumaver etal.i 1201(11 iRusli etaLll2010l) . Regardless of 
the use of longslit or IFU spectrographs, one crucial issue in 
BH mass measurements is spatial resolution: this must be small 
enough to spatially resolve the regions where the gravitational 
effects of the BH can be disentangled from those of the host 
galaxy. Even with the advent of Adaptive Optics (AO) assisted 
observations the best spatial resolution achievable are of the or- 
der of ~ 0.1" which corresponds to ~ 10 pc at a distance of 20 
Mpc. 

This paper is the second in a series dealing with gas kine- 
matical BH mass measurements based on a new method. This 
method, based on spectroastrometry, provides a simple but ac- 
curate way to estimate BH masses and partly overcomes the 
limitations due to spatial resolution which plague the "classi- 
cal" gas (or stellar) kinematical methods, ei ther using longslit or 
IFU spectra. In the first paper of the series (Gnerucci et al.l20To[ 
hereafter Paper I) we illustrated how the technique of spectroas- 
trometry can be used to measure black hole masses focusing on 
the basis of the spectroastrometric approach and showing with 
an extended and detailed set of simulations its capabilities and 
limits. While we mostly focussed on the application of spec- 
troastrometry to longslit spectra, we also showed the technique 
can be extended to integral field spectra. 

In this paper we apply the spectroastrometric method devel- 
oped in Paper I to estimate the BH mass using real data. As a 
benchmark for our spectroastrometric approach to the study of 
local BHs, we selected the galaxy Centaurus A (NGC 5128) be- 
cause it has been extensively studied with the gas kinematical 
method showing that the gas is circularly rotating and that BH 
mass and other free parameters are well constrained from the 
observed kinematics. Moreover both longslit and IFU data are 
available and this allows a direct comparison of the application 
of spectroastrometry to different kinds of data. 

In Sect. 2 we recap the existing measurements of BH mass 
for Centaurus A. In Sect. 3 we briefly resume the results of Paper 
I on the application of spectroastrometry to rotating gas disks for 
the detection of the central BH. In Sect. 4 we apply the method 
to the longslit ISAAC spectra of the nucleus of Centaurus A. In 
Sect. 5 we apply the method to integral field SINFONI spectra of 
the nucleus of Centaurus A, both with and without the assistance 
of Adaptive Optics. Finally, in Sect. 6 we compare and discuss 
the results from the different datasets, drawing some conclusions 
on the reliability and accuracy of the method. 



2. Previous measurements of the black hole mass 
in Centaurus A 

Existing measurements of the black hole mass in Centaurus A 
are summarized in Fig. [1] where the various measurements are 
shown as a function of the spatial resolution of the observations 
(Full Width at Half Maximum - FWHM - of the Point Spread 
Function - PSF) used for each measurement. In that Figure we 
also show the measurement obtained in this paper as a function 
of the angular resolution actually obtained with the spectroas- 
trometric technique (i.e. a fraction of the PSF FWHM as will be 
discussed in detail in the following). 

The supermassive black hole in Centaurus A was first de- 
tected and its mass measured with a near infrared gas kinemat- 
ical study using seeing limited spect ra obtained with ISAAC 
at the ESO VLT (iMarconi et al.ll2001h . Subsequent higher spa- 
tial resolution gas kinematical studies based on longslit spec- 
troscopy were performed using STIS on the HST (Marconi et al] 
l2006h and AO as sisted observations wit h NAOS-CONICA 
at the ESO VLT iHiiring-Neumaver et alJ (12006.) . More re- 
cent studi es based on int e gral fi eld spectroscopy were per- 
formed by Kraino vic et al.l (12007 ') using seeing limited obser- 
vations with CIRPASS at the Gemini South telescope and by 
Neumaver et al. (2007) using AO-assisted observation s obtained 
with SINFONI at VLT On the other hand. lSilge eTal] (i2005) and 
Cap pellari et al.l (l2009l) performed near infrared stellar kinemati- 
cal studies based, respectively, on seeing limited longslit spectra 
(GNIRS at Gemini South) and AO-assisted integral field spectra 
(SINFONI at the ESO VLT). 

The top panel of Fig. [T]shows the different Mbh values ob- 
tained by the previous authors, spread over almost an order of 
magnitude. To understand the origin of these differences, in the 
bottom panel of Fig. [T]we plot the MsHsin^i values, i.e. the val- 
ues constrained by the observed velocity fields and not depen- 
dent on the inclination of the rotating gas disks. In the case of 
the stellar kinematical studies, the authors assumed edge on ax- 
isymmetri c potentials, therefore no correction is made to obtain 
MBHsin^i: 'Sike et all (l2005l) and ICappellari et"an (12009 ) dis- 
cuss the systematic uncertainties associated with their assump- 
tions of edge on models. After removing the inclination effect all 
gas kinematical measurements show statistical fluctuation within 
two times the respective sigma; as noted several times, the incli- 
nation of the rotating disk is an important source of uncertainty 
in gas kinematical measurements. 



3. The spectroastrometric measurement of black 
holes masses 



In Paper I (iGnerucci et alJl20Toh we illustrated how the tech- 
nique of spectroastrometry can be used to measure the black 
hole masses at the center of galaxies. In that paper we focused 
in explaining the basis of the spectroastrometric approach and in 
showing with an extended and detailed set of simulations how 
this method is able to probe the principal dynamical parameters 

of a nuclear gas disk. 

The spectroastrometrical method (see ' Bailevll 199 8') consists 
in measuring the photocenter of emission lines in different wave- 
length or velocity channels. It has been used by several authors 
to study pre-main sequence binaries and the presence of inflows, 
outflows or the disk structure of the gas surrounding pre-main se- 



quence stars (Takami et al. 2003; Baineset 



oundmg p 
al.l2004 li 



Porter et aP 



2004, 2005, Whelan et al. 2005). We compared fliis technique 
with the standard method for gas kinematical studies based on 
the gas rotation curve and showed that the two methods have 
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Fig. 1. BH mass measurements for Centaurus A from the 
works mentioned in the text (top panel) and the correspond- 
ing MsHsin^i values (bottom panel) as a function of the spatial 
resolution of the observations of each measure ment (see text). 
Note t hat the uncertainties on the measurement bv lMarconi et al.l 
( 1200 lb are reduced because they were including uncertainties on 
i. 



complementary approaches to the analysis of spectral data (i.e. 
the former measures mean positions for given spectrum veloc- 
ity channels while the latter measures mean velocities for given 
slit position channels). The principal limit of the rotation curves 
method resides in the ability to spatially resolve the region where 
the gravitational potential of the BH dominates with respect to 
the contribution of the stars. In Paper I we showed that the fun- 



damental advantage of spectroastrometry is its ability to provide 
information on the galaxy gravitational potential at scales sig- 
nificantly smaller than the spatial resolution of the observations 
(~ 1/10, and better as we will also show in the present work). 

The general principle of the spectroastrometric method and 
its capability in overcoming the spatial resolution limit is illus- 
trated by the following simple example: consider two point-like 
sources located at a distance significantly smaller than the spatial 
resolution of the observations; these sources will be seen as spa- 
tially unresolved with their relative distance not measurable from 
a conventional image. However, if spectral features, such as ab- 
sorption or emission lines at different wavelengths, are present in 
the spectra of the two sources, the light spatial profiles extracted 
at these wavelengths will show the two sources separately. From 
the difference in centroid positions at these two wavelengths one 
can estimate the separation between the two sources even if this 
is much smaller that the spatial resolution. This "overcoming" 
of the spatial resolution limit is made possible by the "spectral" 
separation of the two sources. 

For clarity, we summarize here the principal features and the 
main steps of the method presented and discussed extensively in 
Paper I. 

- From the longslit spectrum of a continuum-subtracted emis- 
sion line one constructs the "spectroastrometric curve" by 
measuring the line centroid along the slit for all wavelength 
channels. The "spectroastrometric curve" of the line is given 
by the position centroids as a function of wavelength. 

- From the simulations presented in Paper I, we showed that 
the information about the BH gravitational field is predom- 
inantly encoded in the "high velocity" (hereafter HV) range 
of the spectroastrometric curve which comprises the points 
in the red and blue wings of the line. The HV part of the 
line spectrum originates from the gas moving at high veloc- 
ities closer to the BH, in the case of Keplerian rotation and 
this emission is spatially unresolved. For this reason the HV 
part of the line spectrum is not strongly influenced by the 
spatial resolution or other instrumental effects like slit losses 
or by the intrinsic line flux distribution. On the other hand 
"low velocity" emission (hereafter LV) is usually spatially 
resolved (i.e. the gas moving at lower velocities is located 
farther away from the BH) and this makes the spectroastro- 



3 



A. Gnerucci et al.: Spectroastrometry of rotating gas disks: II. Application to Centaurus A. 



metric curve not useful. The HV range of the curve is iden- 
tified by measuring the spatial extent of the line emission as 
function of wavelength and discarding the central (i.e. LV 
range) bins where the line emission becomes broader than 
the instrumental spatial resolution. 

- By measuring the spectroastrometric curves of a given line 
from at least three spectra taken at different slit position an- 
gles, one can obtain a "spectroastrometric map" of the source 
on the plane of the sky by geometrically combining the three 
curves. In the case of integral field spectra, this step is obvi- 
ously not necessary because the spectroastrometric map is 
derived directly from the data cubes. In the case of a ro- 
tating disk with a radially symmetric line flux distribution 
the points of the spectroastrometric map should lie on the 
disk line of nodes. However even for IFU data the effect of 
slit losses or a non-symmetric line flux distribution can per- 
turb the light centroid positions moving them away from the 
disk line of nodes. As shown in Paper I and discussed above, 
these effects become negligible for the HV range of the map 
where the emission is spatially unresolved. In contrast, the 
LV points of the map tend to lie away from the line of nodes 
in a typical "loop" shape. The final spectroastrometric map 
is then obtained by selecting only the HV points. 

- One can then estimate the disk line of nodes by a line fitting 
of the HV range of the spectroastrometric map, project those 
points on the estimated disk line of nodes and obtain the disk 
rotation curve. Finally, one can apply a simple model fitting 
procedure and obtain the parameters determining the gas ro- 
tation curve, and in particular the BH mass. 

4. Longslit spectra: observations and data analysis 

4.1. The data 

We use available near infrared spectra of the nucleus of 
Centaurus A jibtained with ISAAC at the ESO VLT telescope 
(see iMarconi et al.ll2006l for details). Briefly, the spectra were 
obtained with a Of.' 3 wide slit and cover a wavelength range 
of 1.24 - l.SOyum with a resolving power of A/AA - 10500, 
corresponding to a spectral resolution of ~ 1.2A at the cen- 
tral wavelength (/l= 1.274 jum). The spatial scale of the spec- 
tra is 0.147"/pixel along the slit axis and the dispersion is 
0.58A/pixel. The spatial resolution of the spectra is estimated 
as ~ O'.'S (FWHM of the PSF). There are three diff'erent spectra 
characterized by a diff'erent slit position angle: the one we will 
refer to as "PAl" has a po sition angle of 32.5°, "PA2" has -44.5° 
and "PA3" has 83.5° (see lMarconi et al.ll2006l for details). 

4.2. Analysis of the spectra 

In Fig. |2] we display an example of a typical spectrum of 
Centaurus A extracted from one pixel along one of the slits: 
in particular, this corresponds to the "PAl" continuum sub- 
tracted spectrum, extracted at the position of the continuum 
peak. Several gas emission lines can be identified: Pa/3 at ~ 
1.284yum, [Fell] and He I at ~ 1.281|/m, [Fell] at ~ 1.259yL(m 
and [SJXJ at ~ 1.255yum. The lines used for the gas kinematics 
bv lMarconi et al.l (12006i) are the Pa/S and the [Fe II] lines and we 
will also concentrate on these. 

As clearly visible in the figure, both the Pa/3 and [Fe II] lines 
are blended on the blue sides with [FeII]H-HeI and [SIX], re- 
spectively. Potentially this constitutes a problem in measuring 
the spectroastrometric curve because the blended lines will cer- 
tainly affect the position of the light centroid along the slit for a 



Original data longslit spectrum 



185 


_ 




[FeII] + HeI 


Pa;8 


_ 


1 80 










175 


It 






, 


170 




1 






165 











1.280 1.282 1.284 1.286 1.288 

Wavelength (/J.m) 

Synthetic reconstructed longslit spectrum 




1.280 1.282 1.284 1.286 1.288 



Wavelength (/^m) 

Fig. 3. Position velocity diagram of the observed ISAAC spec- 
trum at PAl near to the Pa/3 line. Upper panel: observed spectrum 
of the Pa/? and [Fe II] H- He I complex. Bottom panel: "synthetic" 
reconstruction of the deblended Pa/3 spectrum. The horizontal 
dotted line overplotted on each panel represents the continuum 
peak position. The isophotes denote the same values in both pan- 
els. 



given blue velocity and the position centroid will not be indica- 
tive anymore of the mean position of the gas at a given line of 
sight velocity. 

In order to solve this problem it is necessary to deblend the 
lines under examination and therefore we performed a simulta- 
neous fit of the 4 lines, each with a Gauss-Hermite function (the 
red solid line), at all slit positi ons along the PAl, PA2 and PA3 
slits. Following iMarconi et al. 112006) we assume that all lines, 
hence Pa/3 and [Fell], share the same kinematics: in the simul- 
taneous fit all lines are constrained to the same velocity, velocity 
dispersion and Hermite parameters (/13 and h^), while they can 
have different line fluxes. Fig.|2]shows an example of such a fit. 

To obtain the spectroastrometric curve we considered all fit- 
ted profiles of Pa/S and [Fell] and we reconstructed synthetic 
longslit spectra of the emission lines, cleaned in terms of noise 
and blended lines. In Fig. [3] we show an example of one origi- 
nal Pa/3 longslit spectrum compared with its "synthetic" version 
where one can notice that the noise has been smoothed away 
and the [Fe II]-t-HeI complex has been removed. It should be no- 
ticed that whereas we constrained Pa/) and [Fe II] to the same 
kinematics, their fluxes can be different, therefore, in terms of 
the spectroastrometric analysis, the two lines are distinct and can 
still provide different results. Additionally the "synthetic" recon- 
structed spectrum is noise-free because each row represents the 
fitted parametric profile but one has to take into account the er- 
rors on the free parameters in order to estimate the errors on each 
synthetic pixel counts. Therefore, for each spectral fit, we simu- 
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Fig. 4. Spectroastrometric curves of the Pa/? (black points) and [Fe II] (red points) lines (ISAAC spectra) at the three slit position 
angles. Left panel: PAl. Central panel: PA2. Right panel: PA3. The dashed vertical lines on each panel represent the limits of the 
HV range. 



lated 1000 synthetic spectra from 1000 realization of the set of 
the five parameters distributed following a pentavariate distribu- 
tion with the fit correlation matrix. The flux and error of each 
pixel is then estimated from mean and standard deviation of the 
1000 realizations. We note that for each profile along the slit, 
the eiTors on fluxes are uncoiTelated because they originate from 
independent fits to the line profiles. 

From the "synthetic" Pa/? and [Fe II] spectra for the three slits 
(PAl, PAl and PA3), following the method outlined above, one 
can derive the spectroastrometric curves which are shown in Fig. 
m Wavelengths are converted in velocity using as reference (zero 
velocity) the rest frame Pa/J and [Fe II] wavelengths (respectively 
1.28216jumand 1.25702//m). 

Fig. ID reveals features which needs further comments. The 
extremely large error bars observed in the Pa/? curve at PAl 
{AQQkml s < v < 55Qkm/s) are due to the fact that the light pro- 
file is spatially resolved and in this range our method of measur- 
ing the centroid position is not reliable. As observed in section 
2.3 this is expected in the LV range, but the relevant informa- 
tion are concentrated in the HV range of the curve. Except for 
those points, errors on photocenter positions range from ~ 0.01 
to ~ 0.05" that is from ~ 1 /50 to ~ 1 / 10 of the spatial resolution 
of the data and this is the accuracy with which we can measure 
centroid positions. 

As expected, the spectroastrometric curves for the two lines 
are marginally different. Their differences are due to the intrinsic 
flux distribution of the lines on the sky plane but, as observed 
in section 2.3.4, these differences tend to disappear in the HV 
range. 

4.3. The spectroastrometric map of the source 

For each emission line we have obtained three spectroastromet- 
ric curves, one for each PA of the slit. Each spectroastrometric 
curve provides the photocenter position along one slit, i.e. the 
position of the photocenter projected along the axis defined by 
the slit direction. Combining the spectroastrometric curves we 
can thus obtain the map of photocenter positions on the plane of 
the sky for each velocity bin. In principle, the spectroastrometric 
curves from two non-parallel slits should suffice but we can use 
the redundant information from the three slits to recover the 2D 



sky map as described in Section 4. 1 of Paper I. We have chosen a 
reference frame in the plane of the sky centered on the center of 
PAl slit (that correspond to the position of the continuum peak 
along the slit) with the X axis along the North direction. For a 
given velocity bin we then determined the position of the light 
centroid on the sky plane resulting in the 2D spectroastrometric 
map shown in Fig.|5] 

Note that the coordinates on the plane of the sky of the cen- 
ter of the PA2 and PA3 slits must be considered as free parame- 
ters. These unknowns are estimated simultaneously with the po- 
sition of the photocenter following ax^ minimization procedure. 
The final spectroastrometric map on the plane of the sky is that 
given by the best fitting set of slit centers. The error bars on 
the points represent the uncertainties resulting from the fit. The 
black points correspond to the HV range (i.e. v < 380 km/s and 
V > 800 km/s) which were actually used to determine the loca- 
tion of slit centers. Indeed, in Paper I we concluded that the HV 
range of the spectroastrometric curve is more robust, and less 
affected by slit losses which artificially change the photocenter 
position in the LV range. 

As observed in appendix A of Paper I, we used the width of 
the Gaussian fitted to the principal peak of the light profile to se- 
lect the HV range: in Fig.|6]we display the FWHMs (Full Width 
Half Maxima) of the Gaussian fitted to the principal peak of the 
light profile for spectroastrometric curves of Centaurus A. We 
can observe that in the LV range the FWHM increases because 
the emission peak is spatially resolved. We compare the FWHM 
to the spatial resolution (~ 0.5") because the FWHM of the light 
profile of an unresolved source should be of the order of the spa- 
tial resolution. We selected the HV range by imposing that the 
FWHM is lower than 1 . 1 times the spatial resolution (FWHM of 
the PSF), resulting in v < 380km/i and v > 800km/i for the 
spectroastrometric curves of both lines. 

The two-dimensional spectroastrometric map just described 
(and shown in Fig.|5]l can now be used to estimate some geomet- 
rical parameters of the nuclear gas disk. If the gas kinematics 
are dominated by rotation around a point-like mass (the BH), 
the position of the light centroid in the HV range should lie on 
a straight line (which identifies the direction of the disk line of 
nodes) and should approach, at increasing velocities, the posi- 
tion of the BH. These considerations allows us to estimate the 
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Fig. 6. FWHM of the Gaussian fitted to the principal peak of the 
light profile for the ISAAC emission line spectra. Upper panel: 
Pa/? line. Bottom panel: [Fell] line. Solid lines: PAl curves. 
Dotted lines: PA2 curves. Dot-dashed lines: PA3 curves. The 
horizontal dashed lines denotes 1 . 1 times the PSF FWHM. 



position angle of the line of nodes (Olon) by fitting a straight 
line to the HV points, considering their errors in both X and Y 
directions (see Fig. |5j- The results of these fits are reported in 
Table [1] together with the formal fit uncertainties. We verified 



the reliability of these unc ertainties using the bootstrap method 
(lEfron & Tibshiranill 19941) . In particular, we randomly extracted 
100 datasets from the HV points and re-performed the fitting of 
the line of nodes. Due to the random extraction, the new datasets 
will have the same number of points as the original one but with 
some points replicated a few times and others entirely missing; 
we thus randomly assign different weights in the fits of the HV 
point. After performing the 100 fits, we estimate the error on 
Slon by taking the standard deviation of the best fit values which 
are usually normally distributed. This error is consistent with the 
formal fit error 

We can also make a first estimate of the BH location by tak- 
ing the average position of the HV points in the spectroastromet- 
ric map. These positions are then refined with the model fitting 
procedure described below and are reported in Table[T] Estimates 
from different lines are all consistent with each other when tak- 
ing into account the ~ 0.01" uncertainties ( ~ 1 /20 of the spatial 
resolution of the data). 

In top panels of Fig. |9]we show the derived spectroastromet- 
ric maps for the Pa/3 band [Fell] lines. It should be noticed that 
all LV points lie outside of the line of nodes, as expected. 

4.4. Estimate of the BH mass from the spectroastrometric 
map 

Here we recover the BH mass value from the spectroastrometric 
map, following the method outlined in Sect. 5 of Paper I. 

Briefly, under the assumption that the gas lies in a thin disk 
configuration inclined by / with respect to the plane of the sky 
(/ = face-on) and the disk line of nodes has a position angle 
Olon, the circular velocity of a gas particle with distance r from 
the BH is given by: 



G[MBH + M/L-L(r)] 

Vrot-^l (1) 

where r is the distance to the BH, L{r) is the radial luminosity 
density distribution in the galactic nucleus and M/L is the mass 
to light ratio of the stars (see Paper I for more details). 
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The component of Vi-ot along the Une of sight (hereafter Vd, 
for "channel velocity") is: 



Vch = Vror Sin(/) + Vsy.s 



(2) 



where / is the inclination of the disk an d we als o added the sys- 
temic velocity of the galaxy V,,., (see iMarconi et aL,2006i, for 
details). 

As explained in section l4!2l and 1431 and sections 3 and 5 of 
Paper I, we only make use of the HV points of the spectroastro- 
metric 2d map to estimate the BH mass. 

The first step is to recover the disk line of nodes by fitting a 
straight line on the 2d map. We then project the position of the 
2d map points (xch,ych) on the line of nodes, calculating their 
coordinate with respect to this reference axis (Sch) and then their 
distance r from the BH used in Eq.[T](i.e. r - k\S d, - So\ where 
So is the coordinate along the line of nodes of the BH and k is 
a scale factor to transform arcsec in the right distance uniQ; see 
Paper I for details). 

Then from Eqs.[T]and|2]we obtain the model channel veloc- 
ity: 



„ ^ G(MBH + M/L-L(k\Sdr-SQ\) . ^.^ 
Vch = ± aI — — s'n(') + Vsxs (3) 



k\S 



ch 



■s. 



and the sign depends on the side of the disk considered. The 
unknown parameters of this model are found minimizing the 
quantity 



X 



Vrh - V, 



ch 



AiSciuPar) 



(4) 



where A{Scii',par) is the uncertainty of the numerator As 
previously discussed, we restrict the fit (i.e. the sum over the ve- 
locity channels) to the HV range. As explained in sect. 5 of Paper 
I, the channel velocity Vdum has no associated uncertainty since 
it is not a measured quantity but the central value of the velocity 
bin. Finally we add a constant error (A,,.,) in quadrature to the 
quantity A{Sch', par) in order to obtain a reduced close to 1 
(see sect. 5 of Paper I for a detailed explanation of this choice). 
Finally, using the best fit values of the model parameters we can 
compute the {x,y) position of the BH in the sky plane. 

We verified that the BH mass value is insensitive to the ac- 
tual value of the position angle of the line of nodes. Indeed, when 
repeating the fit with 9lon values varying within the uncertain- 
ties, Mbh changes by only ~ 0.02dex, well below the typical Icr 
uncertainties of the fit (see Table[T]). We have also checked if the 
final BH mass estimate could be biased by incorrect estimates of 
both 9lon and the BH position. To do this, we repeated the anal- 
ysis including only the HV points which deviate by more than 
2cr from the linear fit to the line of nodes. The final BH mass 
changes by less than 0.1 dex. 

In the rotation curve model shown in Eq. [3] disk inclination 
i and BH mass are coupled since coordinates along the line of 
nodes iSch) do not depend on / and this parameter appears only 
as a scaling factor on the velocity. The reason for this coupling is 
that we are effectively measuring velocities of rotating material 
located on the line of nodes, thus removing any dependence on 
i except for the projection of the velocity along the line of sight. 
Therefore our fitting method can only measure Mbh sin^ /, and 



' To be consistent with all previous measurements we assume a dis- 
tance to Centaurus A of 3.5Mpc. At this distance 1" corresponds to 
~ \lpc. 



we need to assume a value for the inclination to obtain a value 
of the mass. In conclusion the free parameters in our fit are: 

Mbh sin^ / mass of the BH; 
MIL sin^ / mass to light ratio of the nuclear stars; 
line of nodes coordinate of the BH; 



systemic velocity of the galaxy; 



in the following we will only report M sirp-i values and we will 
discuss the inclination values we assume and the relative mass 
values we obtain. 

We have performed the fit of the spectroastrometric data 
from the {Fell} and Pa^S lines, both separately and simultane- 
ously. Fit results are tabulated in Table [T] and presented graphi- 
cally in Fig. |2]where we plot the r = -5o| vs. y = |Vc/i- V.sv.sl 
rotation curve and the line of nodes projected rotation curve (e.g. 
Sch vs. Vch)- The solid red lines represent the curves expected 
from the model (r vs. \Vch - Vjyj| and Sch vs. Vch respectively). 

The first result is that assuming a disk inclination of 25° as in 
IMarconi et aP (1200 6) we obtain Mbh = 10'^ '^^" " ^MfTi, perfectly 
consi stent with the result presented in that work (IMarconi et al.l 
12006' report Mbh = 10^ ''^='° °'*Mfr,). Since we are actually us- 
ing the same ISAAC data of IMarconi et al.l 12006. we can draw 
the important conclusion that the spectroastrometric method pro- 
vides results perfectly consistent with the rotation curve method, 
but with a much simpler approach which does not require to take 
into account complex instrumental effects. 

Furthermore, we note that MjL is not constrained, a re- 
sult which has been already found with the classical rotation 
curve fitting. Computing the radius of the BH sphere of influ- 
ence for Centaurus A and using a stellar velocity dispersion of 
CTstar - I50km/s (Giiltekin et al. 2009) one obtains rBH - 25pc, 
a factor of 20 larger than the distances we are considering here 
(at the distance of 3.5Mpc 1" corresponds to ~ 11 pc and so 
the apparent dimension of rBH is ~ 1.5")- This implies that at 
these small radii the contribution of the stellar mass to the grav- 
itational potential is negligible and consequently M/L cannot be 
constrained with the fit. 

Another important result derived from the application of the 
spectroastrometric method is that the minimum distance from 
the BH at which there is a velocity estimate is ~ 0.05" corre- 
sponding to ~ 0.8 pc, while with the standard rotation curve 
method the minimum distance from the BH which can be ob- 
served is of the order of half the spatial resolution (0.25"). This 
clearly shows how spectroastrometry can overcome the spatial 
resolution limit. In Centaurus A the sphere of influence of the 
BH is already well resolved in ground based observations with 
good seeing and this enables us to obtain a BH measurement 
already from the standard methods. The potential of spectroas- 
trometry clearly relies on the possibility of using the extra res- 
olution to measure the masses of BHs with smaller spheres of 
influence, like those with lower masses, or located in more dis- 
tant galaxies. 



5. Integral field spectra: observations and data 
analysis 

5.1. Data and spectra analysis 

We use available near infrared spectra of the nucleus of 
Centaurus A obtained with SINFONI at the ESO VLT (see 
Neumayer et al. 2007 for details). In particular we make use of 
H band spectra observed in seeing limited mode and H and K 
band spectra obtained with the assistance of the adaptive optics 
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Fig. 8. Example of the FWHM of the fitted 2D gaussian to each 
velocity channel map for the H band seeing limited SINFONI 
[Fell] data. The horizontal dashed lines denotes 1.1 times the 
PSF FWHM. 



(AO) system MACAO. The seeing of the observations as mea- 
sured by the seeing monitor was FWHMy ~ 0.5" (transformed 
to K band , FWHMk ~ 0.38"). Seeing limited spectra use a 
pixel scale of 0.125" x 0.250" and cover a 8" x 8" field of view. 
The SINFONI AO module used as guide star an R ~ 14 mag 
star 36" southwest of the nucleus providing a spatial resolution 
of ~ 0.12" (FWHM). AO assisted spectra use a pixel scale of 
0.050" X 0.100" and cover a 3" x 3" field of view. Spectral res- 
olution is R ~ 4000 for the K band and a slightly lower for the 
H band, R ~ 3000. The total on-source exposure for the K band 
data cube was 13500i whereas for H band was 3600i. For all 
d etails of observations a nd data reduction the read er is referred 
to lNeumaver et al.ll2007l and lCappellari et al.ll2009i 

All data cubes were continuum subtracted by fitting a power 
law function to the spectrum of each spatial pixel (with emis- 
sion and absorption lines masked) which was then subtracted. 
As observed in sections 3 and 4. 1 of Paper I, for this applica- 
tion of spectroastrometry it is mandatory to use continuum sub- 
tracted spectra. The focus is exclusively on gas kinematics and 
the presence of an underlying continuum can significantly alter 
the spectroastrometric measurement of the emission line gas by 
modifying the spatial light distribution in each velocity channel. 
In the following spectroastrometric analysis we will use [Fell] 
observed in the H band (1.6468yL(m) and the H2 line observed in 
K band (2.1259yL(m). Wavelengths are then converted in velocity 
using as reference (zero velocity) the rest frame [Fe II] and H2 
wavelengths (respectively 1.6435yum and 2.1213jUot). 

5.2. The spectroastrometric map of the source 

As observed in section 6 of Paper I, the extension of the spec- 
troastrometric technique to integral field spectra is straightfor- 
ward, and deriving a 2d spectroastrometric map becomes triv- 
ial due to the 2d spatial coverage of integral field spectra. The 
analysis of data now reduces to fitting a 2d Gaussian to each 
channel map in turn, yielding the X, Y positions of the photo- 
centers as a function of velocity, i.e. the spectroastrometric map. 
Therefore, we can directly derive the 2d spectroastrometric map 
from the continuum subtracted SINFONI data cubes, overcom- 
ing the problems related to the uncertainties in slit positioning 
with respect to the galaxy nucleus. 

To select the HV range in the case of integral field data we 
also make use of the widths of the 2d Gaussians fitted to each 



velocity channel map. As explained in sect. 3 and appendix B 
of Paper I we search for unresolved spatial emission. In Fig. IS] 
we show as an example the FWHMs of the fitted 2d Gaussiaqj 
for the seeing limited [Fell] data in the H band where the esti- 
mated seeing is ~ 0.4". We can observe an evident peak in the 
FWHM in the LV range due to the presence of spatially resolved 
emission. As previously, we identify the HV range by consider- 
ing FWHMs lower than 1.1 times the spatial resolution (~ 0.4"). 
The resulting range in this particular case is v < 250km/. s and 
V > 950km/5. 

In the lower panels of Fig. |9]we show the derived spectroas- 
trometric maps for the H band [Fe II] data (both AO assisted and 
not) and for the K band H2 AO assisted data. Typical uncertain- 
ties in the light centroid positions are of the order of ~ 0.01" for 
the seeing limited data (spatial resolution ~ 0.4" and pixel scale 
0.125" X 0.250"), and a factor ~ 2.5 lower for the AO assisted 
data (~ 0.004" with spatial resolution of ~ 0.12" and pixel scale 
of 0.05" X 0.10"). It should also be noticed that all LV points lie 
outside the line of nodes, as expected. 

As previously explained, these 2d spectroastrometric maps 
can now be used to obtain the direction of the disk line of nodes 
and a first estimate of the BH position on the plane of the sky 
(see Table [T]i. The derived BH positions are consistent within 
~ 0.02". The origin of the sky plane corresponds to the contin- 
uum peak position measured directly on the SINFONI datacubes 
with uncertainties of ~ 0.006" for the seeing limited data and 
~ 0.002" for the AO data. Taking into account these system- 
atic errors in setting the origin of the map and the typical uncer- 
tainties of the BH position estimate, BH positions are consistent 
among each other. 

We also note that the derived 0lon values for the [Fe II] line 
observed with or without the AO are consistent within ~ 8°. The 
differences between AO and seeing limited observations might 
also originate because of a warping in the gas disk which results 
in different ori entations of the line of nodes at different spatial 
scales (see also lNeumaver et al]|2007h . 

From Fig.|9]we can also observe that the BH position derived 
from ISAAC J band spectra are shifted by ~ 0.08" with respect 
to that derived from SINFONI H band spectra. This effect can 
be accounted from the fact that the continuum peak positions 
measured in J and H band (i.e. the maps coordinates origin) can 
be different due to the different effect of dust reddening on the 
two bands. 

5.3. Estimate of the BH mass from the spectroastrometric 
map 

Here we estimate the BH mass using the spectroastrometric map 
as described in Sect. 5 of Paper I. The application is exactly the 
same as in the case of the longslit spectra presented in section 
14.41 because in both cases we use as input the 2d spectroastro- 
metric map (i.e. the positions of the light centroids in the plane of 
sky Xci„ych as a function of the corresponding channel velocity 

Vch). 

We have performed the fit of the data from the three spec- 
troastrometric maps of Fig. |9] Fit results are tabulated in Table 
[T]and presented graphically in Fig. [TOl where we plot either the 
r - \S ch - S o\ vs. V = \Vch - Vsysl rotation curve and the line of 
nodes projected rotation curve (e.g. Sch vs. Vd,). The solid red 



- Actually we fit a non circularly symmetric 2D Gaussian function. 
In Fig. [8] we show the minimum of the two FWHM values along the 
proper axis. 
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Table 1. Fit Results. 



Parameter 






Best fit value±error 








ISAAC data 




SINFONI data 




[Fell] line fit 


Pap line fit 


both lines fit 


NO AO [Fell] line fit 


AO [Fell] line fit 


ewN 


-20.3 ± 1.2 


-18.3 ± 1.8 


-19.8 ± 1.0 


-47 ± 2 


-39.2 ± 0.8 


logioiMsH sin-i/MQ) 
logio(M/L sin-i/MQ) 

So ["] 
Vj„ [km/s] 
Ajj., [km/ s] 


7.39 ± 0.02 
-10.1 ±0.0 <2> 
0.105 ±0.008 

593.5 ± 15.6 
11 


7.39 ± 0.03 
-9.6 ± 0.0 <2) 
0.103 ±0.009 
584.8 ± 17.4 



7.39 ± 0.02 
-7.9 ± 0.0 
0.104 ±0.006 

591.4 ±9.7 



7.51 ±0.06 
-14.5 ±0.0 <2) 
0.054 ± 0.006 
599 ± 33 



7.53 ±0.01 
-9.1 ±0.0 
0.018 ±0.003 
592 ± 19 
19 


Xl, (X^ID.O.F.) 


1.01 (9.11/9) 


0.82 (8.20/10) 


0.93 (21.49/23) 


0.55 (6.59/12) 


1.02 (12.2/12) 


Xbh ["] 
yBH ["] 


0.012 ± 0.006 
0.099 ± 0.008 


-0.003 ± 0.007 
0.098 ± 0.008 


0.007 ± 0.005 
0.098 ± 0.005 


-0.001 ± 0.007 
0.037 ± 0.005 


-0.002 ± 0.002 
0.014 ± 0.002 



Notes. Best fit parameter estimated in the fit of the line of nodes. Parameter not constrained from the fit. 



lines represent the model rotation curves (r vs. \Vch - Vjvsl and 
Sch vs. Vch)- 

From the [Fell] data we estimate a value of the BH mass of 
\ogyQ{MBHsirP'i/MQ) - 7.5, perfectly consistent between AO- 
assisted and seeing limited (tab.[T]) whereas from the H2 line we 
obtain a value (~ 0.2 dex) larger Note that the higher accuracy of 
the light centroid positions of the spectroastrometric map for the 
AO assisted [Fell] data with respect to the seeing limited ones 
results in a lower uncertainty on the MBHsin^i best fit value. 

M/L\s not constrained from the fit even with SINFONI data. 
As observed in section 14.41 the radius of the BH sphere of in- 
fluence for Centaurus A is kbh - 14.9/7C coiTesponding to an 
apparent dimension of ~ 0.9"; here we are studying the rota- 
tion curve at ~ 1/20 smaller scales where the contribution of the 
stellar mass to the gravitational potential is negligible. 

An impressive result of the spectroastrometric method is that 
the minimum radii at which we can probe the rotation curve are 
~ 25mas for seeing limited data and ~ 2Qmas for AO assisted 
data. The latter are only slightly smaller but have a much bet- 
ter positional accuracy, as shown in Fig. \W\ These values are 
~ 1/16 and ~ 1/6, respectively, of the spatial resolution (~ 0.4" 
for the seeing limited and ~ 0.12" for the AO assisted observa- 
tions) and correspond to distances from the BH of ~ QAlpc and 
~ 0.35pc, respectively, ~ 1 /40 of the radius of the BH sphere of 
influence. This is a clear demonstration of the great potentials of 
spectroastrometry in overcoming the spatial resolution limit. 

5.4. SINFONI H2 line spectra 

As described above, we also have available AO assisted 
SINFONI spectra in K band where the molecular hydrogen line 
(H2) is observed. Neumaver et al. (2007) used this emission line 
to study the nuclear gas kinematics and to estimate the BH mass. 
They find that this line spectrum has a good S/N and spatial res- 
olution and that the derived gas kinematics are well explained 
by a rotating disk. In contrast the kinematics from the [Fe II] H 
band spectra clearly show disturbances from the presence of the 
jet. 

However, H2 proved to be of lower quality for the purpose 
of spectroastrometric analysis. In particular, the line emission 
is spatially well extended and the S/N in the HV line wings is 
much lower with respect to the [Fell] line. At first sight this 
mi ght seem surprising giv en the high quality of the data obtained 
by Neumaver et al.ll2007l but this is due to the different features 



of the spectroastrometric method. In the spectroastrometric ap- 
plication, as extensively discussed in Paper I, we search for a 
spatially unresolved emission in the HV wings of the line as a 
"signature" of the gas kinematics of the BH gravitational po- 
tential and discard all the spatially resolved emission in the LV 
range. On the contrary, in standard gas dynamical studies, it is 
important to concentrate on the spatially resolved LV range of 
the line spectrum, and particular gas kinematics features like the 
presence of inflows, outflows or jets can be detected in this con- 
ditions. 

In Fig.[TT]we show the comparison between the line profiles 
of the [Fe II] and H2 lines. Both spectra extracted from the re- 
spective SINFONI datacubes (AO assisted observations) from a 
circular aperture of ~ 0.75" centered on the continuum peak. 
Clearly the H2 line spectrum has almost no signal in the HV 
range compared to the [Fe II] line. For this reason the application 
of our method results in a poor quality spectroastrometric map. 
On the other hand in the present work we find the SINFONI 
[Fe II] spectra very useful and it seems our results are not dom- 
inated by the contamination of the jet as observed for the [Fell] 
kinematics by Neumaver et aL2007> This is because we are con- 
centrating on the unresolved HV gas emission that comes from 
the inner region of the nuclear disk and that seems not to be 
strongly influenced by the presence of the jet. 

Summarizing, [Fell] is detected both on large spatial scales 
(those probed with the classical kinematical analysis) and on 
small spatial scales (those probed with the spectroastrometric 
analysis). On large spatial scales and low velocities, the [Fell] 
kinematics is likely affected b y the presence of the j et, while 
this is not the case for H2 (see iNeumaver et al]|2007l for more 
details). On small spatial scales and high velocities, the [Fell] 
kinematics is little affected by the presence of the jet, if any, 
while H2 emission is absent and therefore cannot be used for 
spectroastrometry. This different behavior of [Fell] emitting gas 
is might be related to the warping of the disk, which at small 
spatial scales tends to be perpendicular to the jet axis. 

6. Discussion and conclusions 

We have obtained new mass measurements for the nuclear BH 
of the Centaurus A galaxy by applying the spectroastrometric 
method to longslit and both seeing limited and AO coiTected 
IFU spectra. Compared to the standard method based on rotation 
curves analysis, the spectroastrometric method is much simpler 
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Table 2. Mass estimates for Centaurus A 



"Classical" method applications 


logioiMsH 




ISAAC J band TFe III and Pan lines ("Marconi et al. 2006) (/ = 25°) 


7.39 H 


:0.04 


NACO H band AO TFelll line (Harins-Neumaver et al. 2006) (/ = 45°) 


7.48 


+0.04 


CIRPASS J band Pa^ line (Krainovic et al. 200T) (/ = 25°) 


7.2 


-0.3 


SINFONI K band AO H, line (Neumaver et al. 2007) ((' = 34°) 


7.1 H 


:0.1 


Spectroastrometric method applications 


ISAAC J band [Fe II] line 


7.39 H 


:0.02 


ISAAC J band Pa^ line 


7.39 H 


:0.03 


ISAAC I band [Fe II] and Pa^ simultaneous fit 


7.39 H 


:0.02 


SINFONI H band no AO [Fe II] line 


7.51 H 


:0.06 


SINFONI H band AO [Fell] line 


7.52 H 


:0.01 



15 
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Fig. 11. Comparison between the [Fe II] and H2 line spectra ex- 
tracted from a circular aperture of ~ 0.75" centered on the 
continuum peak (from SINFONI AO assisted observations). 
Overplotted the HV range limits used for the [Fe II] analysis. The 
Okm/s velocity coiTespond to the respective central line wave- 
lengths and the counts are rescaled to approximately match the 
two line peaks. 



from the modeling point of view; it only requires the determina- 
tion of the "spectroastrometric map" and is relatively insensitive 
by the problems that plague the standard approach based on ro- 
tation curves like, e.g., the effect of beam smearing, the intrinsic 
flux distribution of the line, and the biases due to the slit posi- 
tioning in longslit observations. 

With our proposed spectroastrometric approach we can de- 
rive two-dimensional plane-of-the-sky spectroastrometric maps 
of the source characterized by accuracies of position measure- 
ments much lower than the spatial resolution of the observations. 
The mean accuracies of our light centroid position estimates are 
~ 0.01" for the ISAAC data and seeing limited SINFONI data 
(~ 1/40 of the spatial resolution) and ~ 0.004" for AO assisted 
SINFONI data (~ 1/30 of the spatial resolution). The position 
angles of the disk line of nodes estimated from those maps are 

19° for the two ISAAC maps (consistent within ~ 2°). The 

estimates from SINFONI maps are ~ 20° offset from the previ- 
ous. However the SINFONI [Fell] observations (seeing limited 
and AO assisted) provides line of nodes PA estimates consistent 
within -8°. 

The differences between the various estimates of the disk line 
of nodes PA might be due to the different ways in which the 
spectroastrometric maps are constructed, i.e. combining three 
long-slit spectroastrometric curves for ISAAC data and directly 
from the datacube for SINFONI data. Moreover there are many 
clear indications that the Centaurus A nuclear disk is warped 



(seel Neumaver et aill2007t) and therefore the different spectroas- 
trometric maps might probe the average disk line of nodes at 
different spatial scales. In any case, the derived rotation curves 
and hence the BH mass estimates are consistent and not affected 
by the differences in PA. 

One possible cause of concern could be the supp ort provided 
agains t gravity by turbulent motions. However iMarconi et al.l 
(l2006h showed that the line widths of ionized hnes and Pa/S 
are consistent with unresolved rotation and therefore there is 
no indication for turbulent pr essure. Moreove r the BH mass de- 
termination based on H2 by iNeumaver et all i2007) which in- 
cli ides turbulent support is similarly consistent with the result 
bv lMarconi et al.l QOOG) (once the different disk inclinations are 
taken into account, see Fig.[TJ indicating that its effect is small. 

The most important result in this work is the demonstration 
of the capability of spectroastrometry to overcome the spatial 
resolution limit and estimate BH masses in a simple and neat 
way. The minimum distance from the BH at which we can probe 
the gas rotation curve is ~ 5Qmas for the ISAAC data (~ 1/10 
of the spatial resolution) and ~ IQmas for the SINFONI data 
(~ 1/15 of the spatial resolution for the seeing limited data and 
~ 1 /6 for the AO assisted data). In the case of Centaurus A, these 
coiTesponds to ~ 1 /30 and ~ 1/50, respectively, of the radius of 
the BH sphere of influence indicating that is is possible to probe 
deep in the BH potential well, where the contribution from the 
mass in stars is negligible. 

In table |2] we compare MsHsin^i values resulting from re- 
cent applications of the classical rotation curve method to var- 
ious data sets of Centaurus A (including those analyzed here) 
with our results from the application of the spectroastrometric 
method. From Table |2] several conclusion can be reached. 

- Our new simple method based on spectroastrometry is in ex- 
cellent agreement with the classical method based on the ro- 
tation curves, at least when comparing the results obtained 
from the same dataset (cf. the ISAAC Marconi et al. 200^ 
data). This is a fundamental indication for the robustness of 
our new method: using the same dataset we can apply in- 
differently the classical and the spectroastrometrical method 
obtaining perfectly consistent results. 

- When applying the two methods to different datasets but with 
the same target l ine we obtain consistent estimates (cf. the 
Mbh estimate bv lHaring-Neumaver et al.ll2006i targeting the 
H band [Fell] line and the on e presented in this paper, and 
that by Krainovic et al.ll2b07l based on the J band Pa/3, with 
our own which is within ~ 0.2 dex). 

- In general all the measurements reported in Table |2] are con- 
sistent within +0.2 dex and these differences are mainly due 
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to the different data sets and target lines and not from the 
application of a particular method. 

- The application of the spectroastrometric method to differ- 
ent data type (IFU and longslit) give consistent result (within 
only ~ 0.1 dex) and this demonstrates the versatility of our 
method. The application to IFU data is obviously much sim- 
pler but this agreement also show that our method of "recon- 
struction" of the 2d map from multiple longslit spectroas- 
trometry is correct. 

- The application of the method to IFU data with and with- 
out AO also produces consistent results at similar spatial 
scales. This clearly demonstrates how spectroastrometry is 
much less sensitive to spatial resolution than the classical 
method. As expected, the accuracy of measurement posi- 
tions without AO is worse. This is due to the fact that for 
each velocity channel we recover the centroid position by 
fitting a two dimensional Gaussian function and the width of 
this function clearly decreases with increasing spatial resolu- 
tion of the data. Since this function is centro-symmetric we 
can always recover the correct center position. Increasing the 
width of the function only makes the center position uncer- 
tainty larger but does not alter its position. 



A typical feature of spectroastrometry is its insensitivity to 
disk inclination. In our fit, in fact, mass and disk inclination 
are coupled. This is because coordinates along the line of nodes 
(Sch) do not depend on /, which only appears as an unavoid- 
able projecting factor on the velocity. Therefore our method can 
return only a Mbh sirP'i value with the need of assuming an in- 
clination value to derive the value of the mass. In this paper we 
decided to present only the Mbh sirP'i values and to move the 
discussion on the Mbh values to a discussion on the / values 
estimated or assumed by the various authors. For completeness 
we report in Table |2] the inclinations estimated or assumed in 
those previous works. Summarizing, we have applied the spec- 
troastrometric method presented in Paper I to several datasets of 
the nucleus of the Centaurus A galaxy obtaining Mbh estimates 
which are consistent with the classical method based on rotation 
curves. 

In conclusion, the application to Centaurus A has shown that 
the spectroastrometric method is much simpler and straightfor- 
ward than the classical rotation curves method and provides the 
major advantage of enabling spatial scale much smaller than the 
spatial resolution to be probed. It is thus possible to obtain more 
robust and accurate BH mass measurements where the classi- 
cal method fails. Therefore the spectroastrometric method can 
be used to extend the measured Mbh range to smaller BHs in the 
local universe and to similar BHs in more distant galaxies than is 
currently possible. Moreover this method can be applied to any 
type of long-slit or integral field spectra without any particular 
regard for instrument or wavelength range. A very interesting 
and promising application, for example, could be in the sub-mm 
high spatial resolution data provided by ALMA. 

So far the spectroastrometric method has been applied inde- 
pendently from the classical method. However as discussed in 
Paper I, it is clear that the spectroastrometric and "classical" ro- 
tation curves are complementary and orthogonal descriptions of 
the position velocity diagram. Therefore, a future development 
of this method will be its application in combination to the clas- 
sical method based on rotation curves. This will also allow us to 
constrain the disk inclination and thus remove mass-inclination 
degeneracy. 
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Fig. 7. Results of the spectroastrometric modeling of the ISAAC [Fe 11] data (upper panels), ISAAC Pa/? data (central panels) and 
simultaneous fit of the {Fell} and Pa/3 data (lower panels). The left panels show the line of nodes projected rotation curve Sch vs. 
ychan- The right panels show the r - \Sch - Sol vs. V - \Vcii - Vsysl rotation curve. The solid red lines represent the curves expected 
from the model 
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Fig. 9. 2D spectroastrometric map for the ISAAC longslit J band data (upper panels): Pa/3 line (upper left panel) and [Fe II] line (up- 
per right panel). Spectroastrometric map for the SINFONI H band data (lower panels): [Fe II] line for the seeing limited observation 
(lower left panel) and AO assisted [Fe II] line (lower right panel). The red point marks the inferred BH position. The red solid line 
represents the line of nodes of the disk obtained with a linear fit of the HV points, as described in the text. The dotted red lines 
represent the Icr uncertainties on the line of nodes position angle. All boxes have the same angular dimension on the plane of sky 
(~ 0.3" X 0.3"). 
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Fig. 10. Results of the fit for the SINFONI seeing hmited [Fe II] data (upper panels) and SINFONI AO assisted [Fe II] data (lower 
panels). The left panels show the line of nodes projected rotation curve Sd, vs. Vch- The horizontal and vertical dotted lines denote 
respectively the BH line of nodes coordinate and the systemic velocity Vsys- The right panels show the r = \Seh - So\ vs. 
V = \Vch - Vsysl rotation curve. The solid red lines represent the curves expected from the model. 
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